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Abstract 

The study of fractional variational problems with derivatives in the 
sense of Caputo is a recent subject, the main results being Agrawal's nec- 
essary optimality conditions of Euler-Lagrange and respective transver- 
sality conditions. Using Agrawal's Euler-Lagrange equation and the 
Lagrange multiplier technique, we obtain here a Noether-like theorem 
for fractional optimal control problems in the sense of Caputo. 
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1 Introduction 

Noether's theorem, published in 1918 [23], is a central result of the calculus of 
variations that explains all physical laws based upon the action principle. It 
is a very general result, asserting that "to every variational symmetry of the 
problem there corresponds a conservation law" . Noether's principle gives pow- 
erful insights from the various transformations that make a system invariant. 
For instance, in mechanics the invariance of a physical system with respect 
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to spatial translation gives conservation of linear momentum; invariance with 
respect to rotation gives conservation of angular momentum; and invariance 
with respect to time translation gives conservation of energy [201, EI] • The cal- 
culus of variations is now part of a more vast discipline, called optimal control 
[25J , and Noether's principle still holds in this more general setting [30~1 132] . 

Fractional derivatives play an increasing role in mathematics, physics and 
engineering [U [161 El ED EH EH] and the theory of the calculus of variations 
has been extended in order to deal with more general systems containing non- 
integer derivatives [H [21 El El El Ell ESI ESI EZ] • The new fractional variational 
calculus provide a more realistic approach to physics [S1E1E1ES1E21ES], permit- 
ting to consider nonconservative systems in a natural way — a very important 
issue since closed systems do not exist: forces that do not store energy, so- 
called nonconservative or dissipative forces, are always present in real systems. 
Nonconservative forces remove energy from the systems and, as a consequence, 
the standard conservation laws cease to be valid. However, it is still possible to 
obtain a Noether-type theorem which covers both conservative and nonconser- 
vative cases [TOlEJ]. In a more general way, formulations of Noether's theorem 
were proved for fractional problems with left and right Riemann-Liouville frac- 
tional derivatives E31 El ES] . 

In [3] Agrawal proves a version of the Euler-Lagrange equations for frac- 
tional problems of the calculus of variations in the sense of Caputo. One of 
the interesting aspects of the new theory is that both Caputo and Riemann- 
Liouville derivatives play a role in Agrawal's Euler-Lagrange equations. Here 
we use the results of [3J to formulate a Noether-type theorem in the general 
context of the fractional optimal control in the sense of Caputo. 

2 Fractional Derivatives 

In this section we collect the well-known definitions of fractional derivatives in 
the sense of Riemann-Liouville and Caputo (see [H [31 [2H EU EH])- 

Definition 2.1 (Fractional derivative in the sense of Riemann-Liouville). 

Let f be an integrable continuous function in the interval [a, b] . For t £ 
[a, b], the left Riemann-Liouville fractional derivative a Dff{t) and the right 
Riemann-Liouville fractional derivative t D^f(t), of order a, are defined by 




(1) 




where n £ ~N, n — 1 < a < n, and T is the Euler gamma function. 
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Definition 2.2 (Fractional derivative in the sense of Caputo). Let f be an 
integrable continuous function in [a,b]. For t 6 [a,b], the left Caputo frac- 
tional derivative ^Dff(t) and the right Caputo fractional derivative fD£f(t), 
of order a, are defined in the following way: 

'•WV-T^fr-^iff '<■<>», (3) 
= r<sbj jy ~ ^ H) " mM • (4) 

where n E N, n — 1 < a < n. 

Remark 2.3. I/aGN, equalities (JTJ) — (SJ) give the classical derivatives 

aD?f(t)= C a D?f{t)= (!) "/(*), 

«)=?W=(-^°/(t)- 

Remark 2.4. T7ie Caputo fractional derivative of a constant is always equal 
to zero. This is not the case with the Riemann-Liouville fractional derivative. 



3 Main Results 

Our main result is a Noether-type theorem for fractional optimal control 
problems in the sense of Caputo (Theorem 13. 1 7j) . As a corollary, we obtain 
a Noether theorem for the fractional problems of the calculus of variations 
(Corollary DD. 

The fractional optimal control problem in the sense of Caputo is introduced, 
without loss of generality, in Lagrange form: 

I[q(-),u(')]= [ L{t,q{t),u{t))dt— »min, (P c ) 

J a 

c a D?q{t) = ip{t A {t) ) u{t)) , 

where functions q : [a, b] — > M. n satisfy appropriate boundary conditions. The 
Lagrangian L : [a, b] x W 1 x R m — > R and the velocity vector ip : [a, 6] x R" x 
W 71 —>■ W n are assumed to be functions of class C l with respect to all their 
arguments. We also assume, without loss of generality, that < a < 1. In 
conformity with the calculus of variations, we are considering that the control 
functions u(-) take values on an open set of M m . Along the work we denote by 
diL, i = 1,2,3, the partial derivative of function L(-,-,-) with respect to its 
ith argument. 
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Definition 3.1 (Process). An admissible pair (g(-) , u(-)) which satisfies the 
control system ^D^q{t) = tp(t,q(t),u(t)) of problem <\Pq\) is said to be a pro- 
cess. 

Remark 3.2. Choosing a = 1, Problema (\Pc\) is reduced to the classical 
problem of optimal control theory [25]: 

I[q(-),u(-)} = j L(t,q(t),u(t))dt^ mm, (5) 

J a 

q{t) = ip{t,q{t),u{t)) . 

Remark 3.3. The fundamental fractional problem of the calculus of varia- 
tions in the sense of Caputo, first introduced in |3]/ 

I[q(-)) = f L(t,q{t), c a D?q{t))^mm, (6) 

J a 



is a particular case of (\Pc\) ■' we just need to choose cp(t, q, u) = u. 

The fractional Hamiltonian formalism introduced in [11] is easily adapted 
to our present context. Using the standard Lagrange multiplier technique, we 
rewrite problem (\Pp\j in the following equivalent form: 

I[q(-),u(-),p(-)}= I [H(t,q(t),u(t),p(t))-p(t)-°D?q(t)]dt^mm, 

J a 

(7) 

where the Hamiltonian H, is defined by 

H (t, q, u,p) = L (t, q,u)+p-(p (t, q, u) . (8) 

Remark 3.4. In the context of classical mechanics, p is interpreted as the 
generalized momentum. In the optimal control literature, the multiplier p is 
known as the adjoint variable. 

We now proceed with the usual steps for obtaining necessary optimality 
conditions in the calculus of variations. We begin by computing the variation 
51 of functional (J7j): 

51= f [d 2 n ■ 5 q + d 3 n ■ 5u + d 4 n -Sp-Sp- c a D« q - P -s ( c a D« q )] dt , (9) 

J a 

where 5q, 5u and 5p are the variations of q, u, and p respectively. Using the 
results obtained in [HE1E], equation ([9]) is equivalent to 

51= f [{d 2 H - t D%p) ■ 5q + d 3 H ■ 5u + (d 4 H - °D?q) ■ 5p] dt 

J a 



Fractional optimal control and Noether's theorem 



5 



where the fractional derivatives of p(t) are in the sense of Riemann-Liouville 
(in contrast with the fractional derivative of q(t) which is taken in the sense 
of Caputo). Standard arguments conduce us to the following result. 

Theorem 3.5. If (q(-),u(-)) is an optimal process for problem Q-PcD ? then 
there exists a function p(-) G C l ([a, b]; M. n ) such that for all t e [a, b] the tuple 
(<?(■), u(-),p(-)) satisfy the following conditions: 

• the Hamiltonian system 

d 2 H(t,q(t),u(t),p(t))= t D?p(t), 
d 4 H(t,q(t),u(t),p(t))= °D?q(t); 

• the stationary condition 

d 3 H{t,q{t),u{t),p{t)) = 0; 

• the transversality condition 

with Tt given by (JED ■ 

Definition 3.6. A triple (q(-), u(-),p(-)) satisfying Theorem \3.5\ will be 
called a fractional Pontryagin extremal. 

Remark 3.7. For the fundamental fractional problem of the calculus of 
variations in the sense of Caputo ([6]) we have Ft = L + p ■ u. It follows from 
Theorem 1 3. 51 that 

t D%p = d 2 L, (10) 
d 3 H = & p = -d 3 L t D b a p = - t D£d 3 L . 

Comparing both expressions for t D^p, we arrive to the fractional Euler- Lagrange 
equations proved by 0. P. Agrawal §<?/: 

d 2 L+ t D%d 3 L = 0. (11) 

In other words, for the problem of the calculus of variations ([6]) the fractional 
Pontryagin extremals give Agrawal 's Euler-Lagrange extremals. 
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Remark 3.8. Our optimal control problem Q-PcD only involves Caputo frac- 
tional derivatives but both Caputo and Riemann-Liouville fractional derivatives 
appear in the necessary optimality condition given by Theorem \3.5[ This is dif- 
ferent from FJJ [7^1 where the necessary conditions only involve the same 
type of derivatives (Riemann-Liouville) as those in the definition of the frac- 
tional optimal control problem. This fact also occurs in the particular case of 
the calculus of variations as noted in J3]/: the Riemann-Liouville derivative is 
present in the necessary condition of optimality (TTTj) but not in the formulation 
of the problem ([6]). 



The notion of variational invariance for problem QPcD is defined with the 
help of the equivalent problem ([7]). 



Definition 3.9 (Invariance of Q-PcP without transformation of time). We 



say that functional (\Pq\) is invariant under the one-parameter family of in- 
finitesimal transformations 

q(t) = q(t) + e£~(t, q, u,p) + o(e) , 
u(t) = u(t) + eq(t, q, u,p) + o(e) , 
p{t) = pit) + eg(t, q, u,p) + o(e) , 



if and only if 



[H{t,q{t),u(t),p(t))-p(t) ■ C a D«q{t)] dt 

[H(t,q(t),u(t),PV)) ~P(t) ■ aD?q(t)] dt (12) 



for any subinterval [t a ,tb] Q [a, b]. 

Noether's theorem will be proved following similar steps as those in [T3J for 
the problems of the calculus of variations in the Riemann-Liouville sense. 

Lemma 3.10 (necessary and sufficient condition of invariance). If func- 
tional (\Pc\l is invariant, in the sense of Definition \3.9[ then 



d 2 H (t, q, u, pH+W (t, q, u, p).q+ (d,H (t, q, u, p) - C a D?q)-o-p- c a D?t = . 

(13) 

Proof. Since condition (fl2l) is to be valid for any subinterval [t a , h] C [a, b), we 
can write ffl2l) in the following equivalent form: 

H(t,q(t),u(t),p(t)) -p{t) ■ C a D« q {t)=H{t,q{t),u{t),p{t))-p{t) ■ G a D<?q{t) ■ 

(14) 
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We differentiate both sides of ( 1T41) with respect to e and then substitute e by 
zero. The definition and properties of the Caputo fractional derivative permit 
us to write that 

= d 2 H {t, q,u,p)-Z + d 3 H (t, q,u,p)-<; + {d 4 H (t, q, u, p) - ^D^q) ■ g 
d 



n-a-l ( d 

J e=0 



^ de 



= d 2 U (t, q, u, p) ■ f (i, q) + d 3 H (t, q,u,p)-q+ (d 4 H (t, q, u, p) - ^D^q) ■ g 

-"■f^jfc-^ 1 (!)"«*■ (15) 

Expression (TT5|) is equivalent to (|T3|) . □ 

In this work we propose the new notion of fractional conservation law in 
the sense of Caputo. For that, we introduce the operator T>f. 

Definition 3.11. Given two functions f and g of class C 1 in the interval 
[a,b], we introduce the following operator: 

vnf,g] = -gtDtf + f c a D^ g , 

where t G [a, b] and u G Rq". 

Remark 3.12. Similar operators were used in /IT]/ and [TR Def 19]. While 
here T>£ depends both on Caputo and Riemann-Liouville derivatives, the ones 
introduced in 07] [73]/ involve Riemann-Liouville derivatives only. We note that 
in the classical context uj = 1, and V][f, g] = gf + f g' = -| (fg) = V\[g, f). 

Definition 3.13 (fractional conservation law in the sense of Caputo - cf. 
Def. 23 of [13]). A quantity Cf (t,q(t),^D"q(t),u(t),p(t)^ is said to be a frac- 
tional conservation law in the sense of Caputo if it is possible to write Cf as 
a sum of products, 

C f (t,q(t),°D?q(t),u(t),p(t)) 

r 

= £)C* (t,q(t)CD?q(t),u(t),p(t)) ■ Cf (t,q(t)CD?q(t),u(t),p(t)) (16) 
i=i 

for some r G N, and for every i = 1, . . . , r the pair C\ and Cf satisfy one of 
the following relations: 

V? [Cl (t,q(t)CD?q(t),u(t),p(t)),Cf (t,q(t)°D?q(t),u(t),p(t))] = (17) 
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or 

V? [Cf {t,q(t)°D?q(t)Xt),p(t)),Cl (t,q(t)°D?q(t)Mt),P(t))] = (18) 

along all the fractional Pontryagin extremals ( Definition 1 3. 6\) . 

Remark 3.14. If a = 1 f|T7|) and ffT8l) coincide, and Cf f|T6|) satisfy the 
classical definition of conservation law: 4 [Cf(t, q, q)} = 0. 

Lemma 3.15 (Noether's theorem without transformation of time). If func- 
tional < \Pq\ i is invariant in the sense of Definition ^. % then 

is a fractional conservation law in the sense of Caputo. 

Proof. We use the conditions of Theorem 13.51 in the necessary and sufficient 
condition of invariance (TIBl) : 

= - dzH-Z-W-S- (d.H-^D^q) -g + p. C a D«£ 

□ 

Definition 3.16 (Invariance of Q-PcD )- Functional (\Pq\) is said to be invari- 
ant under the one-parameter infinitesimal transformations 



't = t + er(t, g, u, p) + o(e) , 
q(t) = q(t) + e£(t, q, u,p) + o{e) , 
u(t) = u(t) + e$(t, q, u,p) + o(e) , 

j(t) = pit) + eg(t, q, u,p) + o(e) , 



(19) 



if and only if 



[H {t, q{t),u(t),p(t)) - p(t) ■ C a D?q(t)] dt 
-f(tb) 



t(t a ) 



[H (i, q(t),u(i),p(i)) - pit) ■ C a D?q{i)] dt (20) 



for any subinterval [t a ,tb] Q [a, b]. 

Next result provides an extension of Noether's theorem 
optimal control problems in the sense of Caputo. 



for fractional 
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Theorem 3.17 (Noether's theorem for fractional optimal control prob- 
lems). If the functional Q-PcD is invariant under the one-parameter infinitesi- 
mal transformations (fT9l) . then 

C f (t,q(t)°D?q(t),u(t) lP (t)) 

= [H(t, q (t),u(t), P (t))-(l-a)p(t)- c a D?q(t)]T-p(t)-Z (21) 

is a fractional conservation law in the sense of Caputo (see Definition \3.13\) . 

Remark 3.18. If a = 1 problem ( [PcD takes the classical form ([5]) and 
Theorem \3.17\ gives the conservation law of fSDj : 

C{t, q{t), u(t),p(t)) = [H(t, q(t),u(t),p(t))} r - p(t) ■ £ . 

Proof. Every non-autonomous problem flPcD is equivalent to an autonomous 
one by artificially considering t as a dependent variable. For that we consider 
a Lipschitzian transformation 

[a, b]3t i — > o-f(X) E [<7 a , a b ] 

satisfying the condition t' a = = /(A) = 1 for A = 0, such that ([7]) takes 
the form 

I[t(.),q(t(.))Mt(-)),p(t(-))] = / [H(t(a),q(t(a)),u(t(a)),p(t(a)) 

-p{t{a))- c a D<? {a) q{t{a))\ tja , 

where t(a a ) = a, t(at) = b and 

*D?(*)q(t(o-)) 

1 r°fW { rl \ n 



Tin — a) J a \ds 

(4) 

(Cr a C x D:q(a), ( X 



(C 

Then, we have 

I[t(.),q(t(.)),u(t(.)),p(t(.))] 
f ab r 

= / H(t(a),q(t(a)),u(t(a)lp(t(a))-p(t(a)) ■ (t a )~<* ^D a a q(a 

J da 

= [ 1 H f (t(a),t a ,Q(t((T)),u(t(a)),p(t(a)^D:q(a)) da 

b 

[H(t,q(t),u(t),p(t))-p(t) ■ C a D?q{t)] dt 
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If functional I[q(-),u(-),p(-)} is invariant in the sense of Definition 13.161 then 
the functional I[t(-),q(t(-)),u(t(-)),p(t(-))] is invariant in the sense of Defini- 
tion [33J Applying Lemma [3.151 we obtain that 



C f [t(a),C q(t(a)),u(t(a)),p(t(a)),^q(a)) = p(t(a)) • £ + ^{t{a))r (22) 
is a fractional conservation law in the sense of Caputo. For A = 0, 

p(t(a))=p(t) (23) 

and it follows from the stationary condition of Theorem 13.51 (see third equality 
in (HDD) that 



_ &Hf_ _d_ 

^~ dt' ~ dt' 



(t CT ) 2 



-H + p^D?q 

= - ap(t(a)) ■ j\ (a - S )— 1 (£) " q(s)ds -H + p- c a D?q 

(t a ) 2 

= -{H-(l-a)p- c a D?q) . 

(24) 

Substituting (1231) and (1241) into (1221) . we obtain the fractional conservation law 

m- □ 



As a corollary, we obtain the analogous to the main result proved in pT3| 
for fractional problems of the calculus of variations in the Riemann-Liouville 
sense. 

Definition 3.19 (variational invariance for ©)• Functional ([6]) is said to 
be invariant under the one-parameter family of infinitesimal transformations 

i=t + er(t,q) + o(e), 
q(t) = q(t)+eC(t,q) + o(e), 



if and only if 

T L (t, q(t), c a Dfq{t)) dt = f_ {tb) L (i, q(t)°D?q(t), ) dt 

Jt a Jt{t a ) 

for any subinterval [t a ,tb] C [a, b\. 
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Corollary 3.20 (Noether's theorem for fractional problems of the calculus 
of variations). If functional (EI) is invariant under the family of transformations 
(1251) . then 

C, (t, g, c a D<? g) = d 3 L (t, q, c a D«q) • £ 

+ [L (t, g, a c Z^) - ac^L (t, «, C a D?q) ■ c a D<?q] r (26) 

is a fractional conservation law in the sense of Caputo. 

Proof. The fractional conservation law (I26p is obtained applying Theorem 13. 171 
to functional ([SD- □ 

Remark 3.21. If a = 1 problem (jBJ) is reduced to the classical problem of 
the calculus of variations, 

I[q(-)}= [ L(t,g(t),g(*))^min, (27) 

J a 

and one obtains from Corollary \3.20\ the standard Noether's theorem [23]: 

C{t, q, q) = d 3 L (t, g, g) ■ £(t, q) + [L (t, q, q) - d 3 L (t, q, q) ■ g] r(t, q) (28) 

is a conservation law, i. e. (1281) is constant along all the solutions of the Euler- 
Lagrange equations 

d 2 L(t,q,q) = j t d 3 L(t,q,q) (29) 
(these classical equations are obtained from (11 ip putting a = l). 



4 Illustrative Examples 

In classical mechanics, when problem (|27|) does not depend explicitly on g, i.e. 
L = L(t,q), it follows from (fTUl) and (|2"9~1) that the generalized momentum p is a 
conservation law. This is also an immediate consequence of Noether's theorem 
[23] : from the invariance with respect to translations on q (r = 0, £ = 1), 
it follows from (128]) that p = d 3 L is a conservation law. Another famous 
example of application of Noether's theorem in classical mechanics is given by 
the conservation of energy: when the Lagrangian L in (12"7|) is autonomous, i.e. 
L = L(q,q), we have invariance under time-translations (r = 1, £ = 0) and it 
follows from (jSJ), (]T0]) and (128]) that the Hamiltonian TC (which is interpreted 
as being the energy in classical mechanics) is a conservation law. Surprisingly 
enough, we show next, as an immediate consequence of our Theorem 13 .171 that 
for the problem ( \Pq\) with a fractional order of differentiation a (a ^ 1), the 
following happens: 
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similarly to classical mechanics, the generalized momentum p is a frac- 
tional conservation law when L and tp do not depend explicitly on q 
(Example ED); 



differently from classical mechanics, the Hamiltonian Ti is not a fractional 
conservation law when L and (p are autonomous (Example 14.2ft . 



In situation (ii), we obtain from our Theorem 13 .171 a new fractional conservation 
law that involves not only the Hamiltonian Ti but also the fractional order of 
differentiation a, the generalized momentum p, and the Caputo derivative of 
the state trajectory q (see ( 130]) below). This is in agreement with the claim 
that the fractional calculus of variations provide a very good formalism to 
model nonconservative mechanics [Til 127] . In the classical case we have a = 1 
and the new obtained fractional conservation law (I3"0l reduces to the expected 
"conservation of energy" TC. 



Example 4.1. Let us consider problem fl-PcD with L(t,q,u) = L(t,u), 



(p(t,q,u) = tp(t,u). Such a problem is invariant under translations on the 
variable q, i.e. condition (120]) is verified for i = t, q(t) = q(t) + e, u{t) = 
u(t) and pit) = pit): we have dt = dt and condition f )20l) is satisfied since 

?fl?«(*) = ?fl?9(t)-- 



Fiji — a) 
c ;D?q{t) + c a D?s 



= C a D?q(t)- 

According with f|T9|) one has £ = 1 and r = <j = g = 0. It follows from 
Theorem \3.11 that p(t) is a fractional conservation law in the sense of Caputo. 



Example 4.2. We now consider the autonomous problem (\Pq\) '■ L(t, q, u) = 
L(q,u) and tp{t,q,u) = tp(q,u). This problem is invariant under time transla- 
tion, i.e. the invariance condition (j2DJ) is verified for i = t + e, q(t) = q(t), 
u(t) = u(t) and pit) = pit): we have dt = dt and (|20|) follows from the fact 
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1kat°D?ffl = °D?q(t): 





1 




r 


n — 


a) 




1 




r 


n — 


a ) 




1 




r 


n — 


a) 



" t+£ ( d \ n 

it + £ - ey-*- 1 (-) q(9)de 

( d \ n 

(*"*) ( ^) Q(t + e)ds 
= c a D?q{t + e) = C a D?q{i) 
= C a D?q(t) ■ 

With the notation ( |T9i) one has r = 1 and £ = c; = g = 0. We conclude from 
Theorem \3.1T\ that 

H(t,q,u,p)-(l-a)p-CD?q (30) 



is a fractional conservation law in the sense of Caputo. For a = 1 (1301) repre- 
sents the "conservation of the total energy": 

H(t,q(t),u(t),p(t)) = constant , t 6 [a, b] , 

for any Pontryagin extremal (q(-) , u(-) , p(-)) of the problem. 
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